Energy distributions of field emitted electrons from carbon nanosheets: manifestation 

of the quantum size effect 
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We emphasize the importance of experiments with vohage dependent field emission energy distri- 
bution analysis in carbon nanosheets. Our analysis shows the crucial influence of the band structure 
on the energy distribution of field emitted electrons in few-layer graphene. In addition to the main 
peak we found characteristic sub-peaks in the energy distribution. Their positions strongly depend 
on the number of layers and the inter-layer interaction. The discovery of these peaks in field emission 
experiments from carbon nanosheets would be a clear manifestation of the quantum size effect in 
these new materials. 

PACS numbers: 79.70.-|-q, 81.05.Uw, 73.43.Cd 



Recently, freestanding carbon nanosheets (CNSs) have 
been synthesized on a variety of substrates by radio fre- 
quency plasma enhanced chemical vapor deposition [l|, Q ■ 
The sheets are consisting of several graphene layers and 
stand roughly vertical to the substrate. It has been 
found that CNSs have good field emission characteristics 
with promising applications in vacuum microelectronic 
devices [1, 0, [s,, 3] ■ High emission total current at low 
threshold field enables using CNSs as an effective cold 
cathode material. 

Until now only the current-voltage characterization 
was used in studies of CNSs. At the same time, voltage 
dependent field emission energy distribution (V-FEED) 
analysis is known as a powerful experimental method to 
interrogate the field emission. As compared to classical 
I-V characterization, V-FEED analysis can provide more 
information related to both inherent properties of the 
emitter and to the basic tunneling process T\ . In particu- 
lar, in single- walled carbon nanotubes (CNTs) the FEED 
has shown characteristic peaks originated from the sta- 
tionary waves in the cylindrical part of the nanotube [sj . 
Their number and sharpness were found to increase with 
the length of the tubes. Notice that short periodic varia- 
tions were also observed in the thickness-dependent field 
emission current from ultrathin metal films (UMF) 
The calculated electron energy distribution curve char- 
acteristic of UMF was found to have "steps" which cor- 
respond with the quantized "normal" energies [l^. The 
resonant-tunneling peaks with specific microscopic tun- 
neling mechanisms were also observed in field emission 
from nanostructurcd semiconductor cathodes [ll| . A dif- 
ferent example of the quantum size effect in CNTs, which 
originates from the intrinsic properties of the energy band 
structure, was revealed in field emission [1^. It is rea- 
sonable to expect manifestation of quantum size effects 
in subnanometer CNSs. 

In this Letter, we calculate the FEED of electrons from 
CNSs. For this purpose, we take into account the en- 



FIG. 1: The location of a graphene sheet with respect to the 
electric field. 



ergy band structure of few- layer carbon systems resulting 
from the tight-binding approach. Both the field emission 
current (FEC) and the FEED are calculated by using 
the independent channel method suggested recently in 
Ref. [l^. Our analysis clearly shows that the FEC only 
measurements give incomplete information. We found 
that the FEED enables determination of the number of 
layers in few-layer graphene as well as direct verification 
of the high sensitivity of the band structure to the num- 
ber of layers in few-layer graphene reported recently in 
Ref. [11. 

Let us consider the graphene layer in the presence of 
the external electric field F directed along the z-axis (see 
Fig.l). The emitted current density takes the following 
form: 
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where e is the electric charge, h = 2TTh the Planck con- 
stant, e the energy, p momentum, /(e) = [exp(£/fcT) + 
1]~^ the Fermi-Dirac distribution function, D{e,Px,Py) 
the transmission probability of an electron through a po- 
tential barrier, and Vg — de/dpz the group velocity. The 
integrals are over the first Brillouin zone with account 
taken of the positivity of Vg . 

For a two-dimensional (2D) structure, one can use 
the relation J f{px)dpx = f{0)h/lx- Moreover, when a 
graphene sheet has the finite size in the y-direction, py is 
quantized. Therefore, the current density in Eq. ([T]) can 
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FIG. 2: The Brillouin zone of graphite. 



be written as 



2e 

hlxly 



E 



(2) 



where the Fermi energy is chosen to be zero. Limits s'}^^^ 
and e^j„ come from the expUcit form of the band struc- 
ture. 

We suggest that the transmission probabiUty is given 
by the WKB approximation in the form 



D{e) = cxp 



F 



(3) 



where C = 87r(2m)i/2/3eft,, y = (eF/47reo)^/^/0, </> is the 
work function, Eq the dielectric constant, m the electron 
mass. The function v{y) describes a deviation of the 
barrier from the triangle form due to image effects and 
can be approximated as v{y) « 1 — y^-^^ (see Ref. [isj). 

The band structure of graphene multilayers has been 
obtained within the tight-binding approach in Ref. p^ . 
Besides, an approximation to the dispersion relation 
can be found from Slonzewski-Weiss-McClure (SWMcC) 
model for graphite with Bernal stacking SWMcC 
model describes the wave-vector dependence of electron 
energy in the vicinity of the HKH edge of the Brillouin 
zone (see Fig. 2). The electron energy spectrum is ob- 
tained from the equation 



detliJ-el = 0, 
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and 



H = 



Ei = A + ^iT+hiT^, 
^2 = A-7ir+|75r2, 

E3 = i72r', 

^^13 = 7f (-70 + 74r) exp(iQ;)cr, 
H23 = 73(70 + 74r) exp(ia)cr, 
-^33 = 73r exp(ia)cr 



(4) 



(5) 



(6) 



with r = 2cos(fcj_c), a = k\\y/3/2a = p\\Vf/^i, kj_ being 
the wavevector projection onto the direction HKH, fc|| 



the modulus of the wavevector in the j/z-plain, a the an- 
gle between k|| and the direction TK, c the distance be- 
tween nearest neighbour layers, a the lattice constant, P|| 
the momentum in the yz-plain, and vp the Fermi veloc- 
ity. Parameters 7^ describe interactions between atoms 
and A is the energy difference between two sublattices 
in each graphene layer. In Ref. graphite parame- 
ters were estimated as 70 = 3.16 6^,71 = 0.39 6^,72 = 
-0.020 eV, 73 = 0.315 eV, 74 = -0.044 eV, 75 = 0.038 eV 
and A = —0.008 eV. The spectrum of few-layer graphene 
can be obtained from Eq. Q by replacing F by 



F" = 2 cos 



N + 1 



1 , 
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where N is the number of layers. For graphene bilayer 
N = 2 and 72 = 75 = so that Eq. ([4]) with account 
taken of Eq. ([7]) gives the result of Ref. 20] while for 
A'^ = 1 (only 70 differs from zero) it reproduces the known 
tight-binding spectrum of graphene. 

As a first approximation one can neglect all interac- 
tions except between the nearest-neighbor atoms in the 
same layer and between A-type atoms between adjacent 
layers (which are on top of each other), i.e. all parame- 
ters except for 70 and 71 are putted to be zero. Then the 
spectra of multilayers can be approximated by 



e" = ± 



(v/(7lV2j^+^-7i72), (8) 



where 7" = 71F". 

The FEC and the FEED {P{e)) arc connected by (see, 
e.g., Ref. 0) 



dePie). 



(9) 



The explicit form of P{e) for few-layer graphene can be 
found from Eq. Indeed, for layers of a large (infi- 

nite) size the sum in Eq. ([2|) can be replaced by the in- 
tegral and, correspondingly, one has to use £min{Py) and 
^maxiVy) instead of e^j„ and e'^ax- o^'' case, these 
Py-dependent functions can be easily calculated from Eq. 
([8]) . Finally, we have to change the order of integration 
in Eq. The resuh is 



P{e) = ^fie)D{e) ^ 0(e„)/kk, (10) 

where e„ = |£| + 7", g — Ae/{h?Nc), and 9(e) is the 
Heaviside step function. 

For graphene monolayer one gets 



=N 



P~{s) = ^f(s)D{s)\e\. 



(11) 



Within the Fowler-Nordheim (FN) approximation (weak 
fields and low temperatures) the FEC is found to be 
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FIG. 3: Reduced FEED for one- to four-layer graphene. All 
SWMcC parameters except for 70 and 71 are putted to be 
zero, F = 4 V/nm. The peak height for A'' = 1 is chosen to 
be unity. 




FIG. 4: Reduced FEED for three- and four-layer graphene. 
The whole set of SWMcC parameters is used, F = 4 V/nm. 



where b = exp {-C(l)^/^v{y)/ F) and d = SC(l)^/H{y)/2F 
with t{y) w 1 -I- 0.1107yi-33 (see Ref. [H). For bilayer, 
the FEED is obtained as 
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(13) 



and, correspondingly. 



where Ki{x) is the MacDonald function. Notice that 
only the first term in Eq. (jlSp is significant at weak 
fields. When the interlayer interaction is weak (7id <C 1) 
Eg. p^ passes into Eq. p^ . For large 71c? one gets 



d-yi 



J = 
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Thus, we obtain a standard FN exponent while the pre- 
exponential factor becomes proportional to F^/-^ instead 
of F^ for the FN theory. Fig. [3] shows the FEED for 
different numbers of graphene layers. 

Let us now take into account all possible interactions. 
For this purpose, we use the whole set of SWMcC param- 
eters and put a = 0. The numerical results are presented 
in Fig. m It should be stressed that for iV = 1 and N = 2 
the calculated FEEDs are not sensitive to other interac- 
tion constants and curves are found to be identical to 
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FIG. 5: Schematic illustration of the method of independent 
channels for few-layer structures. Black points indicate peaks 
of emission channels. 



those shown in Fig. [31 When TV increases little shifts of 
the minima relative to the Fermi energy are obtained. As 
is clearly seen, for iV > 1 FEEDs have characteristic sub- 
peaks. The number of peaks and their positions strongly 
depend on the number of layers and the interaction con- 
stants, first of all, 71. There is a pronounced depression 
in FEED at the Fermi energy which would be typical for 
3D gapless semiconductors. 

Fig. [5] gives a clear illustration of our results. At room 
temperatures, the Fermi-Dirac distribution function re- 
stricts the FEED above the Fermi energy, so that the 
electron emission from the valence band dominates. The 
transmission probability decays exponentially with de- 
creasing e. Therefore, at fixed F one can estimate the 
energy range of emitted electrons as 



Ae 



3C01/2 



(16) 



In accordance with Fig. [5] the number of emitting chan- 
nels m at the energy e is defined as m = [Apy{s)/Spy] 
where the brackets [ ] indicate integer part. Generally, 
P{e) - C{e)f{e)D{e) where C(e) = m/ly is the den- 
sity of emitting channels. In CNSs ly is large enough and 
C{e) — > Apy{e)/h. In our case, Apy{e) ^ at |e| ^ so 
that P{e) (see Fig. [3]). Interestingly that taken into 
account Eq. ([5]) one can easily calculate C(e) and, cor- 
respondingly, P{e) for any A'^ without integrations. The 
shape of the FEED in Figs. [3] and [4] directly depends on 
the density of emission channels. When e riches the top 
of the next branch of the spectrum this branch becomes 
" switched-on" thus resulting in a distinctive point in the 
FEED. Evidently, the closer a position of the branch to 
Ae the less pronounced is an additional peak in P{e). 

An important difference from the emission of single- 
walled carbon nanotubes should be mentioned. The di- 
ameters of CNTs are very small thus resulting in a set 
of discrete channels. For metallic CNTs there also ex- 
ists at least one emitting channel at |e| 0. However, 
as distinct from CNSs the density C(e) tends to a con- 
stant value at £ ^ and, correspondingly, the FEED 
exhibits behavior typical for conventional metallic emit- 
ters without any minimum near the Fermi energy. On 
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FIG. 6: FEED for bilayer at difFerent temperatures. The 
whole set of SWMcC parameters is used, F = 4 V/nm. 
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FIG. 7: FEC for bilayer in Fowler-Nordheim coordinates at 
different temperatures. 

the contrary, in semiconducting CNTs the FEED has a 
characteristic gap at the Fermi energy (cf. Ref. plj). 

Similar arguments are valid for the emission from the 
conduction band where, however, the limiting role of the 
Fermi-Dirac distribution is of decisive importance. When 
the temperature grows, the electrons from the conduc- 
tion band become involved in the emission. As is seen 



in Fig. [6l the regime of the so-called thermal field emis- 
sion occurs at high temperatures of emitters. Notice that 
there is a rather symmetrical behavior of curves in both 
bands, which is valid for all considered few-layer struc- 
tures. Fig.[7]shows an influence of the thermal field emis- 
sion on the behavior of emission current in the FN coor- 
dinates. The curves are found to be practically identical 
for any N. At high temperatures one can see marked de- 
viations from the standard FN plot in the region of weak 
fields. Notice that similar deviations were observed in ex- 
periments with as-received CNSs (see Fig. 9 in Ref. [1]). 

In conclusion, we have calculated the FEED for few- 
layer graphene films and found the presence of character- 
istic sub-peaks originated from involving in the emission 
process additional branches in the energy spectrum of 
layered structures. Since the peak positions are directly 
determined by the number of layers the discovery of such 
peaks in the FEED would be a clear manifestation of the 
quantum size effect. Therefore, the experimental studies 
of the FEED for CNSs are very relevant. Furthermore, 
the FEED analysis gives a new experimental tool to es- 
timate the inter-layer interaction constants (along with 
Raman scattering in Ref. |22J and photoemission meth- 
ods in Ref. 23]) and provides important information on 
the concrete types of emitting CNSs as well as allows 
one to identify the number of layers in emitting CNSs. 
For example, the absence of sub-peaks would indicate 
that the emission occurs from monolayer graphene. In 
addition, the emitter temperature is taken into account 
in the FEED via the Fermi-Dirac distribution function 
and can be determined by the half-width and the rel- 
ative height of shapes in the conduction band. Finally, 
using an approach suggested in Ref. [23 | one can measure 
the resistivity of CNSs. 
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